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Table Al.l. (continued) 



Exposure Average Exposure Average 

mJD Time Intensity mJD Time Intensity 
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Table A 1.2. Barnard's Star Log of Observations 
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36 


14910.5 


10167.077 


115, 


.95 


14968.8 


9777, 


.169 


34.95 


14929.3 


10167.1 


73, 


.95 


14985.5 


9916, 


.557 


116.85 


14927.3 


10198.508 


116, 


.85 


14903.2 


9916, 


.57 


76.95 


14925.7 


10198.52 


76, 


.95 


14921.4 
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Table A 1.2. (continued) 





Exposure 


Average 




Exposure 


Average 


mJD 


Time 


Intensity 


mJD 


Time 


Intensity 


9916.581 


75 


14919.7 


10198.531 


75.9 


14928.4 


9935.324 


117 


14882.6 


10237.311 


118.05 


14914 


9935.336 


77.85 


14917 


10237.324 


80.85 


14942.5 


9935.347 


76.05 


14934.9 


10237.335 


56.73 


14931 


9962.134 


117 


14984.8 
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ABSTRACT 



We have observed Proxima Centauri and Barnard's Star with Hubble Space 
Telescope Fine Guidance Sensor 3. Proxima Centauri exhibits small-amphtude, 
periodic photometric variations. Once several sources of systematic photometric error 
are corrected, we obtain 2 milh-magnitude internal photometric precision. We identify 
two distinct behavior modes over the past four years: higher amplitude, longer period; 
smaller amplitude, shorter period. Within the errors one period {P ~ 83*^) is twice 
the other. Barnard's Star shows very weak evidence for periodicity on a timescale 
of approximately 130 days. If we interpret these periodic phenomena as rotational 
modulation of star spots, we identify three discrete spots on Proxima Cen and possibly 
one spot on Barnard's Star. We find that the disturbances change significantly on time 
scales as short as one rotation period. 
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Subject headings: stars: flare — stars: individual (Proxima Centauri, Barnard's Star) 
— stars: late-type — stars: rotation — stars: spots 



1. Introduction 

We present photometry of Proxima Cen and Barnard's Star, results ancillary to our 
astrometric searches for planetary-mass companions ( [Benedict et al., 1997aD . Our observations 
were obtained with Fine Guidance Sensor 3 (FGS 3), a two-axis, white-light interferometer 
aboard the Hubble Space Telescope (HST). [Bradley et al., 1991 provide an overview of the FGS 3 



instrument and Benedict et al., 1994a[ describe the astrometric capabilities of FGS 3 and typical 



data acquisition strategies. [Benedict et al., 1993[ assessed FGS 3 photometric qualities and 



presented the first evidence for periodic variability of Proxima Cen. This latter result was based 
on 212 days of monitoring. Subsequent data exhibited a period of variation very nearly twice the 
original ( Benedict et al., 1994b[ ). Since that report, we have obtained 14 additional data sets for 



Proxima Cen and 12 new sets for Barnard's Star. The primary value of these observations lies in 
their precision, not in their temporal span or aggregate numbers. We have previously determined 
that a 90 sec observation obtained with FGS 3 has a 1 — cr precision of 0.001 magnitude at V = 11 
( Benedict et al., 1993[ ), in the absence of systematic errors. 



In this paper we discuss the data sets and assess systematic errors, including background 
contamination and FGS position-dependent photometric response. We also present a revised 
photometric flat field. We then exhibit and analyze light curves for Proxima Cen and Barnard's 
Star. We find weak evidence for periodic variations in the brightness of Barnard's Star. However, 
Proxima Cen exhibits significant periodic photometric variations, with changes in amplitude 
and/or period. We next interpret these variations as rotational modulation of chromospheric 
structure (star spots and/or plages), and conclude with a brief comparison to other determinations 
of the rotation rate of Proxima Cen. Tables |l] and ^ provide aliases and physical parameters for 
our two science targets. 

We use the term 'pickle' to describe the total field of view of the FGS. The instantaneous 
field of view of FGS 3 is a 5 x 5 arcsec square aperture. Figure || shows a finding chart for the 



Barnard's Star reference frame in the FGS 3 pickle as observed on 6 August 1994. Benedict et al. 



1993 contains a finding chart for the Proxima Cen reference frame. 



2. Data Reduction 
2.1. The Data 



All position and brightness measurements from FGS 3 are comprised of series of 0.025 sec 
samples (e.g., 40Hz data rate), of between 20 and 120 sec or 600 sec duration. Each FGS 
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contains four photomultipliers, two for each axis. We sum the output of all four to produce our 
measurement, S, the average count per 0.025 sec sample, obtained during the entire exposure. The 
coverage for both targets suffers from extended gaps, due to HST pointing constraints (described 
Benedict et al., 1993) and other scheduling difficulties. The filter (F583W) has a bandpass 



m 



centered on 583 nm, with 234 nm FWHM. 

For Proxima Cen the data now include 152 shorter exposures secured over 4 years (March 
1992 to October 1997) and 15 longer exposures (July 1995 to July 1996). Each orbit contains 
from two to four exposures. The longest exposure times pertain only to Proxima Cen observations 
obtained within Continuous Viewing Zone (CVZ) orbits. These specially scheduled orbits permit 
~ 90 minutes on field, during which Proxima Cen was not occulted by the Earth. Appendix 1.1 
gives times of observation, exposure times, and average counts, S, for all Proxima Cen photometry. 

Barnard's Star was monitored for three years (February 1993 to April 1996), and observed 
three times during each of 35 orbits. Exposures range between 24 and 123 seconds duration. 
Appendix 1.2 gives times of observation, exposure times, and average counts, S, for all Barnard's 
Star photometry. 



2.2. Background Light 

We first noted that background contamination might be an issue while assessing the use 
of astrometric reference stars for photometric fiat-fielding. These stars are typically far fainter 
than the primary science targets. Using them to fiat field the Proxima Cen and Barnard's Star 
photometry introduced a strong one-year periodicity (and considerable noise, since they are 



fainter stars). This problem was not identified in [Benedict et al., 1993| , since we had access to 
data spanning less than two-thirds of a year. Figure |2| shows S for two faint reference stars in 
the Barnard field plotted against angular distance from the Sun. These stars appear brightest 
when closest to the Sun. Zodiacal light is a source whose brightness depends on the sun-target 
separation. The fitting function in Figure ^ is 

/ = A + i?sin(^), (1) 

chosen to produce a minimum contribution at 9 = 180°. We find A = 137.1 it 0.4 and 
B = —4.1 lb 0.5 counts per 25ms for an average exposure time of 100 sec. At a 60° elongation 
the contamination amounts io V = 22.5 =b 0.3 arcsec"^. The Barnard field is at ecliptic latitude 
(5 = -|-27°. From a tabulation in Allen (1972) we calculate a signal equivalent to V = 22.1 arcsec"^ 
for zodiacal light at 60° elongation and ecliptic latitude (3 = -|-30°. The agreement supports our 
identification of this background source. 

We present in Figure ^ an average S for these two Barnard reference stars plotted as a 
function of time, uncorrected and corrected for background. These data have been fiat fielded 
using the time-dependent response function discussed in section 2.3.2| (equation ^). Note the 
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reduction in the amplitude of the scatter for the corrected photometry. Presuming Zodiacal Light 
as the source, contamination levels are even less for the Proxima Cen observations at ecliptic 
latitude /3 = —44°, introducing a maximum systematic error of 0.0007 magnitude for a 100 sec 
observation. We conclude that the effects of this component of the background are insignificant 
for Proxima Cen and Barnard's Star photometry. 

Should background determination become more important in the future, we note that during 
an intra-orbit observation sequence the PMT are never turned off. Hence, the HST data archive 
contains PMT measurements taken during slews from one star to the next. The astrometric 
reduction pipeline at the Space Telescope Science Institute has been modified to provide these 
background data automatically. 



2.3. Photometric Flat Fielding 

We explore two kinds of flat fielding; position- and time-dependent. We first assess whether 
or not flat-field corrections are necessary, and, if so, determine their functional form. 



2.3.1. Position- dependent Photometric Response 

Having discovered that background variations contaminate the photometry of faint astrometric 
reference stars, we required an alternative source for flat field data. To maintain the astrometric 
calibration of FGS 3, a star field in M35 has been measured roughly once per month for the 
last four years. Whipple et al., 1995| describe this continuing astrometric Long-term Stability 



(LTSTAB) test. The field, on the ecliptic, and, hence, always observed in one of two orientations 
(Fall or Spring) flipped by 180°, contains bright stars for which background contamination is 
negligible. However, an initial application of a time-dependent flat field based on bright M35 stars 
also introduced a strong one-year periodicity. 

The positions of the three M35 LTSTAB stars within FGS 3 are shown in Figure ^ . The 'eye' 
is bordered by the pickle edge at the two nominal rolls for this field. The central circle (diameter 
^ 3'.8) is accessible by the FGS 3 instantaneous aperture for any HST roll. 

Figure ^ (bottom) presents normalized intensities (I = S(t)/Sai,, where Sav is determined 
from the entire run of data) for the three LTSTAB stars as a function of time. The variation 
of each star has first been modeled by a linear trend. The parameters, intercept (/q) and slope 
(/'), are given in Table ^. The resulting residuals (Figure ^, top) have been modeled with a sine 
wave, while constraining P = 365'?25 days. The residuals have a square-wave periodic structure 
because, rather than a range of spacecraft rolls, there are only two orientations. The resulting 
parameters and errors are given in Table ^. In Figure ^ we plot the amplitude of this side-to-side 
variation against radial distance from the pickle center. For the M35 stars, the further the star 
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from the pickle center, the larger the roll-induced variation. Figure ^ includes several other one 
year period amplitudes; a preliminary result for GJ748 {V ~ 11.1, ecliptic lat /3 ~ +23°) always 
observed in the center of the pickle, the Barnard reference star photometry from Figure ^ corrected 
for background, and photometry of the brightest reference star in the Barnard field (Figure |l|, 
star 36, V ~ 11.5). Figure |^ suggests that within the inscribed circle of Figure || (r < 180"), 
position-dependent photometric response variations should be less than 0.002 magnitude. 

We have also identified one high spatial frequency position-dependent flat field component 
for FGS 3. Light curves for two of the Barnard reference stars evidenced sudden decreases in 
brightness with subsequent return to previous levels. The decrease for reference star 34 was 29%; 
for star 36, 17% . Shown in Figure both decreases occurred in the same location within FGS 3, 
very near the -Y edge. The pickle X, Y of the center of this area is (X, Y) = -25, 627. We estimate 
the size of the low-sensitivity region to be ~ 10 x 10 arcsec. Additionally, Proxima Gen reference 
star observations acquired one year prior to the Barnard's Star reference star observations and 
within a few arcsec of this position showed no decrease, providing additional evidence that FGS 3 
is not suitable for wide-field, precise faint star photometry. 

Evidence that the photometric response may vary locally and randomly with time dissuades 
us from mapping a position-dependent flat field over the entire pickle. However, for bright stars 
{V < 11) observed within ~ 20 arcsec of the pickle center (Figure [l|), these identified systematics 
should produce very little effect. All Proxima Cen and Barnard's Star observations were secured 
within 15 arcsec of the pickle center. 

2.3.2. Time- dependent Photometric Response 

Figure ^ indicates that FGS 3 has become less sensitive with age. For all three LTSTAB 
stars the linear trends (/', Table ^ agree within the errors. The apparent 1% drop in sensitivity 
over 1000 days requires confirmation. Figure ^ presents the time varying normalized intensity for 
two other astrometric program stars observed with FGS 3, GJ 623 and GJ 748. Both M dwarfs 
were observed in pickle center. Comparing the I' in Table ^ and Table ^ the rate of decline in 
brightness for GJ 623 and GJ 748 is identical (within the errors) to that seen in the M35 stars. 

A final concern is that the rate of decline of PMT sensitivity might vary with wavelength. The 
M35 stars (stars 547, 500, and 312 in the catalog of Gudworth, 1971) have 0.12 < B -V < 0.49, 
while GJ 623 and GJ 748 have B — V ^ +1.5. There appears to be no dependence on color. 

The weighted average for five stars from three different fields yields 

FF = 1.131 ± 0.006 + (1.30 ± 0.06) x 10~^mJL» (2) 

as the temporal photometric flat field for the pickle center. 

As an additional test of the reality of this sensitivity decrease, we note that the intensity data 
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for the two astrometric reference stars in the Barnard's Star field shown in Figure ^ have been 
flat-fielded with equation ^. Thus, a total of seven stars from four different fields show similar 
brightness trends, adequate evidence for a sensitivity loss in FGS 3. 

2.4. Photometric Calibration 

All magnitudes presented in this paper are provisional, since a final calibration from F585W 
to V is not yet available. If magnitudes are given, they are derived through 

V = -2.5 log(5) + 20.0349 (3) 

with no color term, where S is the average count per 0.025 sec sample, summed over all four 
PMT. No results are based on these provisional calibrated magnitudes. They are provided only as 
a convenience. 

2.5. Summary: Photometric Error and Photometric Precision 

We have identified sky background (Zodiacal Light), within-pickle response variations, 
and time-dependent sensitivity variations as contributing sources of systematic error for our 
photometry. Since our science targets, Proxima Cen and Barnard's Star, are bright, the effect of 
Zodiacal Light is at most 0.001 magnitude. Since we observe these stars only in the pickle center, 
spatially-induced variations are reduced to about 0.001 magnitude, our claimed per-observation 
precision at y ~ 11. A weighted average of the temporal response of five stars in three fields 
provides a very precise flat fleld whose slope error could introduce at most 0.001 magnitude 
systematic error over 1000 days. (Since we are doing only differential photometry we ignore 
the zero point error in the fiat field.) Combining these sources of error yields a per-observation 
precision of 0.002 magnitude. 

3. Photometric Results 

We apply the fiat field (equation Q) to the Appendix 1.1 and Appendix 1.2 S values and 
plot (Proxima Cen, Figure Barnard's Star, Figure [lO| ) the resulting intensities as function of 
modified Julian Date (JD - 2440000). Our coverage in time is not uniform for either target. There 
are extended gaps in coverage, some due to the HST solar constraint (no observations permitted 
closer than ±50° to the Sun). The largest gap, in 1994 for Proxima Cen, was due to an awkward 
transition from Guaranteed Time Observations to Guest Observer status and a hiatus due to 
suspected equipment problems. 
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3.1. Trends in Brightness 

For Proxima Cen the solid line in Figure |9| indicates an overall trend of increasing brightness 
with time. In units of normalized intensity the rate of change of brightness (1.63 it 0.37 x 10~^) 
is similar to that of the adopted flat field (equation ||). For Barnard's Star (Figure ^) the slope 
of the upward trend in units of normalized intensity is +0.91 it 0.18 x 10~^, again, suspiciously 
similar in absolute value to the adopted fiat-field relation (equation ^) . Since seven stars from four 
different fields exhibit the sensitivity decrease discribed by the fiat field, the Proxima and Barnard 
upward trends are unlikely to be a fiat-field artifact. 

A final caveat: Proxima Cen and Barnard's Star are somewhat redder (Tables |l| and ||) than 
GJ 623 and GJ 748. If the sensitivity loss varies with wavelength (e.g., more sensitivity loss for 
blue than for red wavelengths), it would have to be a very steeply dependent function, showing no 
effect at B-V = +1.5. 



3.2. Proxima Cen 

The flat-fielded photometry for each exposure in each orbit appears in Figure |^. The period 
and amplitude variations evident in Figure ^ will be discussed in section 4. Our total time on 
target, obtained by summing the exposure times in Appendix A.l, was 6^'6. Proxima Cen is a 
flare star (V654 Cen) and these data contain exposures 'contaminated' by stellar flares (marked 
Fl - F4 in Figure |^. We identified these events by inspecting the 40Hz photometric data stream 
for each observation. An example of flare contamination (including a light curve) can be found in 
Benedict et al., 1993| , which discusses a slow, relatively faint (AV < —0.10), and multipeaked flare 



on mJD 8906 (Fl in Figure An explosive flare on mJD 9266 {AV ~ ff^6 in one second, F3 in 
Figure ^) produced astrometric changes in Proxima Cen, analyzed in detail by Benedict et al.J 



1997b| . This spectacular event provided the motivation for the subsequent CVZ observations (cz 



in Figure |9|), each permitting 30 minutes of monitoring for flares. The F4 event at mJD 9368 had 
a relatively small amplitude {AV ~ —0.13), but lasted throughout the entire 130^* exposure, hence 
its large effect on the exposure. Walker (1981) predicts a flare with intensity similar to F3 once 
every 31 hours. Thus, while disappointing, it is not surprising that we captured none as bright 
as the F3 event in our additional 2.5 hours of CVZ on-target monitoring. It may be significant 
that we saw any flares at all, since even the small amplitude events have only a 60% chance of 



occurring during our total monitoring duration. We will discuss this further in Section 5.1. 



Individual observations secured within an orbit and not affected by flaring exhibit an internal 
consistency at the 0.002 magnitude level. 
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3.3. Barnard's Star 

The flat-fielded photometry for each of the three Barnard's Star exposures acquired within 
each orbit appears in Figure |l^. Note that the time scale is exactly that used for Figure ^ to 
facilitate comparison. Again, note that those observations secured within an orbit exhibit an 
internal consistency at the 0.002 magnitude level. We find variations within each orbit, but no 
obvious flaring activity in the Barnard's Star results. Possible period and amplitude variations in 



the Barnard's Star data will be discussed in section 4.2.2 



The scatter within each orbit in Figure ^ is somewhat larger than the previously determined 
([Benedict et al, 1993|) 0.001 magnitude measurement precision. In particular we inspected 



the observations on mJD 9935 and 9994 and found only a slight upward slope during the first 
observation on each date. Since the majority of first observations within each orbit are lower, this 
intra-orbit scatter is most likely an instrumental effect, amounting to about 0.001 magnitude. The 
first observation low bias is sometimes seen in the Proxima Cen data (Figure ^). Leaving all first 
observations uncorrected will only slightly enlarge the formal errors for our per-orbit means. 



4. Analysis 

For subsequent analyses of Proxima Cen, we removed the flare contributions by subjecting the 
per-orbit average to a pruning process. All exposures obtained during each orbit are presented in 
Figure |9|. If one exposure differs by more than 2.5 — 0" from the mean for that orbit, it is removed 
and the mean recalculated. This process results in 71 normal points with associated dispersions 
(the standard deviations calculated for 2, 3, or 4 exposures in each orbit) for Proxima Cen. No 
exposures were removed from the Barnard's Star series, since no intra-orbit points (shown in 
Figure ^) violated the 2.5 — a criterion. The resulting per-orbit average S values are presented 
as direct light curves in Figure |l3| (Proxima Cen) and Figure |l6| (Barnard's Star). Forming 
these normal points provides per-orbit photometric precision better than 0.002 magnitude for the 
following analyses. 



4.1. Lomb-Scargle Periodograms 



From Figures |9| and 10 we suspect that there are periodic variations in both the Proxima Cen 
and Barnard photometry. To obtain a preliminary identification of these periodicities we produce 
Lomb-Scargle Periodograms ( Press et al., 199^ ) from the per-orbit normal points presented in 



Figure |13| (Proxima Cen) and Figure |16| (Barnard's Star). 

The most statistically significant period in the Proxima Cen periodogram (Figure 0) is 
at P ~ 83"^, with a false-positive probability less than 0.1%. The very small peak at P '"^^ 42^^ 



indicates the relative strengths of the period derived from the first 212 days ( Benedict et al. 
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1993) relative to the higher ampUtude P~ 83"^ period. The false-positive probability for the period 
derived in that paper from only the first 212 days was ~1 %. Since the periodogram provides 
no results for very short periods, we have some concern that we are undersampling a more rapid 
variation. We can rule out a range of periods 2'^ < P < 20*^ from detailed inspection of clusters of 
data near m JD 8840 (Figure |l^) , where we had a series of closely-spaced (in time) observations 
(see Appendix 1.1). 

The periodogram for Barnard's Star is shown in Figure |l^. The strongest peak(at P ~ 130*^) 
has a 10% false-positive probability. We have much less compelling evidence of variability for 
Barnard's Star than for Proxima Cen. 



4.2. Light Curves 

4-2.1. Proxima Cen Light Curve 

Given strong support for a periodic variation (periodogram, Figure ^l]) and for an overall 
trend in the brightness (Figure we model the per-orbit average variations seen in the direct 
light curve (Figure |l3|) with a sin function and trend 

2-K 

I = I^ + l't + Asm{{—)t + (/)), (4) 



To reconcile the earlier results ( Benedict et al., 1993| ) with the newer data, we first attempted 



to model the entire light curve with only two distinct segments, grouping segments B, C, and 
D together. From the earliest data (segment A) the Proxima Cen photometric variations are 
characterized by a shorter period and smaller amplitude. Later data are best fit with a longer 
period and larger amplitude variation, as evidenced by the periodogram (Figure ^). Parameters 
for these fits are listed A and BCD in Table | (hues 1 and 2). We find Pbcd/Pa = 1-97 ± 0.04. 

Noting very large residuals for segment C, we next explored the possibility that Proxima Cen 
repeats a low to high amplitude cycle by fitting the four segments (A - D, Figure |l^) with the 
same model (equation |^. The parameters for these fits are presented in Table ^ (lines 1 and 3-5). 
Within the errors. Pa = Pc and Pb = Pd, with Pq/Pq = 1.99 ± 0.02. A reduction in the number 
of degrees of freedom by 17% (fitting 71 data points with twenty, rather than ten parameters), 
reduced the residuals by ~ 30%. This relative improvement is some support for alternating high 
and low amplitude states. It is also evident that segments A and C have very nearly half the 
period of segments B and D. 



Figure |14| contains phased light curves for the four segments labeled in Figure 13. In the 
top panel we show that the phase shift required to align the two long-period segments is small 
(A(/» = +0.11) for the period P = 83.5*^ suggested by the periodogram (Figure ^l]). We have 
shifted segment D down by AS = -104.4 counts. The bottom panel of Figure |l^ shows a phased 
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light curve for the two shorter period segments (A and C), phased also to P = 83. 5^^. Shifts 
in phase and intensity to achieve alignment are indicated in the figure. The clean double sin 
wave also demonstrates that the low-amplitude segments, A and C, have half the period of the 
higher-amplitude, longer-period segments, B and D. 

Finally, Figure ^ is used to demonstrate that the same low-amplitude, short-period variations 
seen in segments A and C may also be present in segments B and D. We fit a sin wave to the 
phased B and D light curve in the lower panel of Figure constraining the period to one cycle. 
The top panel of Figure 15 shows the residuals to that fit. These residuals are then fit with a 
sin wave with the period constrained to one-half cycle. Comparing with the bottom panel of 
Figure 14, we find a similar double sin wave, nearly identical phase, and an amplitude (A = 25 ib 8) 
close to that reported for segments A and C in Table ||. 



4-2.2. Barnard's Star Light Curve 

We turn now to the per-orbit average photometry of Barnard's Star. Figure [l^ (bottom) 
contains the per-orbit average direct light curve. The error bars are 1 — a, obtained from the 
dispersion of the three observations within each orbit on each date (Figure |^). Residuals to a 
linear trend are presented in the top panel of Figure |l^. The sin wave fit to these residuals was 
constrained to have the most significant period detected in the Figure ^ periodogram, P = 13rf4. 
Figure ^ contains a light curve for the trend-corrected Barnard's Star photometry of Figure 16 
(top), phased to P = 130^4. The phased light curve is far less clean than for Proxima Cen. The 



periodogram and Figures [1^ and |17| provide only weak evidence for periodic variation, primarily 
due to the poor sampling. 



Discussion of Photometric Results 



Instruments can impress spurious periodicities on data ( [Kristian, 1991 ). It is comforting that 
we find for all segments of either data set that Psam Pprox- 



Stars have local imperfections in their atmospheres (e.g., the Sun, |Zirin 1988 ). Stars other 
than the Sun have been shown to be spotted, photometrically (dwarf M stars; Kron, 1952| ) and 
spectroscopically (e.g., Hatzes, 1993| ; Neff et al., 1995 ). Other M stars have been shown to have 
spots, both dark ( [Bouvier et al., 19951) and bright (a Ori, [Cmiland fc Dupree, 1996|) . 



A spot on a rotating star is a model rich in degrees of freedom. Spots can be bright (plages) 



or dark (see Pettersen et al., 1992 for a discussion of the choice between dark spots on a bright 
background or bright spots on a dark background). Spots can wax and wane in size, driving the 
mean brightness level of a star up or down (IVrba et al.,^988D . Spots can migrate in latitude, 
which, when coupled with presumed differential rotation, can change the phasing and perceived 
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rotation period. Spots are thought to migrate up or down (relative to the star center) within the 
magnetosphere ( Pameron &: Campbell, 1993| ), inducing perceived period changes. In the following 
sections we shall interpret the variations seen in Figures 13 and 16 as rotational modulation of 
spots or plages. 



5.1. Spots on Proxima Cen 

If we assume a fundamental rotation period P = 83. 5"^, then variations in the amplitude 
(Figure |l3|) could be due to spot/plage changes. With the sparse set of single-color photometric 



data presented in Figure 13 we have made no effort to quantitatively model spots (c.f. Neff et al 



19961) . The period and amplitude changes can be qualitatively modeled using plages and spots, 



but require the disappearance of a feature or a major change in feature size or temperature in less 



than one rotation period (e.g., the A to B segment transition seen in Figure 13). 

Segments B and D (Figure |l^, top) require a single large or darker spot that reduces the 
average brightness of Proxima Cen by AV ~ 0.03. This feature was not present in segment A 
and disappeared during segment C. To phase segments B and D we applied a shift of Ai;^ = 0.1 
radians. Thus, the spot site lagged behind the fundamental rotation period by about 5° over the 
end of B to start of D time separation of ~ TOO'^. Whether due to latitude migration coupled with 
differential rotation or from changes of height within the magnetosphere is unknown. If due to 
differential rotation, then either the spot moved very little in latitude or differential rotation on 
Proxima Cen is several orders of magnitude less than the Sun ( ^irin 1988| ). 

Segments A and C exhibit smaller amplitude variations with a period almost exactly half that 



found for segments B and D. These segments (see the phased light curves in Figure |lj, bottom and 
Figure top) could be produced by two smaller spots spaced 180° in longitude, carried around 
by the fundamental 83'?5 rotation period, and persisting through all segments A to D. These two 
spots produce a AV ~ 0.01. One of these small or less dark (warmer) spots lies at nearly the same 
longitude as the prominent spot seen in segments B and D. The other lies near the center of the 
brighter hemisphere in segments B and D. 



From Figure 13 we note that the peaks in segments B and D were brighter. In segment B the 
minima were deeper, implying a darker (cooler) spot. To increase the amplitude of the maxima 
in segments B and D requires the existence of plages, or that the hemisphere not containing the 
single large spot became brighter due to spot changes. One of the pair of spots, associated with 



the brighter hemisphere {(p ~ 0.3, Figure 15, top), does have a shallower minimum than the other 
of the pair at = 0.8. If these are the same spots responsible for the variation in segments A 
and C, then the spot at ~ 0.3 did have a shallower minimum in segments B and D (compare 



Figure |l^, bottom and Figure |l^, top). However, that spot did not become less dark by enough to 
account for the increased maxima seen in segments B and D. As a consequence we propose plage 
activity to increase the segment B and D maxima. 
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Flaring activity seems more prevalent in segment B. As seen in Figure 14, three of four 
flares are grouped near phase (j) = 0.8. Association with this deep minimum might imply some 
connection of flaring activity with the largest or coolest star spot, which is also at the same 
longitude as one of the two smaller spots seen best in segments A and C. The remaining flare, 
F4, lies close to (j) = 0.2, the other spot of the low-amplitude pair of spots. Spot/flare association 
was previously noted in the M dwarf EV Lac, shown also to have longitude-dependent flaring 
associated with a star spot site (|Leto et al., 1997 ). However, for EV Lac, flares were detected 
a year before the spots became easily detectable by their system {/S.V ~ 0.1) and, once spots 
formed, flare activity abated. 



5.2. A Small, Variable Spot on Barnard's Star 

The periodogram (Figure ^) does not provide a clear identification of a single period of 
variation. The trend-corrected direct light curve (Figure |l6[ top) has been fit with a sin wave, 
constraining P = 130?4. The constant amplitude is /S.V ~ 0^01, about five times our formal 
photometric error. Given the sparse coverage, it is speculative to interpret this light curve as 
showing rotational modulation of a single, small spot decreasing in size. 



5.3. Rotation Periods for Proxima Cen and Barnard's Star 



Rotation periods for Proxima Cen have been predicted from chromospheric activity levels by 
Doyle (1987), who obtained P=51"' ± 12^^. Guinan k Morgan (1996) measure a rotation period (P 
= 31^5 =b l'?5) from lUE observations of strong Mg II h+k emission at 280nm. We find no support 
for either rotation period in our periodogram (Figure 11) or light curves, direct (Figures ^) or 
phased (Figure ^4|). We do note that the variation due to the spot pair produces a period between 
the Doyle prediction and Guinan & Morgan measurement. 

The observed variations for Proxima Cen and Barnard's Star, if interpreted as rotationally 



modulated spots, yield rotation periods far longer than for other M stars. For example, [Bouviei 
|et al., 1995 find 4"^ < P < 8^^ for a sample of young, early-M T Tauri stars. Magnetic braking is 



postulated to slow rotation over time ( Cameron &: Campbell, 1993 ). The inferred rotation period 
for Proxima Cen is consistent with old age. This age may be 4 - 4.5 By, if Proxima Cen is coeval 
with a Cen ( Demarque et al., 1986| ). A relatively older age for Barnard's Star can be surmised 
from lower than solar metallicity ( pizis, 1996 ) and higher than solar space velocity (Eggen, 1996), 
both consistent with a longer rotation period, if one accepts the reality of the variation. 
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5.4. Shorter Time-scale Variations 



The level of internal per-orbit precision for these photometric data is near 0.002 magnitudes. 
Hence, the dispersion about the phased light curves for Proxima Cen (Figure 14) is likely intrinsic 
to the stars. Two possibilities are miniflaring and the creation and destruction of small star 
spots and plages. That either phenomenon must have a duration longer than hours, at least 



for Proxima Cen, is suggested by the segment A phased light curve (Figure 14, bottom) and a 



detailed light curve for segment A (see [Benedict et al., 1993| , Fig. 3). Segment A contains four 
pairs of back-to-back orbits and one set of three contiguous orbits (on mJD 8845). In each case 
the time on target coverage is over 90 minutes. For most of these contiguous orbits differences are 
within two standard deviations and not statistically significant. Since 'flare' implies a relatively 
short duration, miniflaring cannot be the cause of the scatter. 



5.5. Activity Cycles 



The ~ llOO'' cycle of alternating high and low amplitude (see Figure 13) is suggestive of 
an activity cycle for Proxima Cen. However, the gap in our coverage in segment C weakens any 
claims that can be made relative to the timing of this cycle. Comparing their 1995 lUE data with 
earlier archival data, Guinan fc Morgan, 1996 propose an activity cycle that was in a low-state in 
1995, agreeing with our identification of segment C representing a low-state (Figure |l^). 



6. Conclusions 

1. For FGS 3 photometry we have identified four sources of systematic error: background 
contamination (primarily Zodiacal Light); spatial flat field variations (significant only for target 
positions r > 20" from pickle center) ; temporal sensitivity changes (calibrated to a level introducing 
a 0.001 magnitude differential run-out error in 1000 days); and a possible warm-up effect (see 
section |3.3| ) . 

2. Two to four short (t < lOO'') observations with FGS 3 during one orbit yield 2 
milli-magnitudes precision photometry, provided the targets are bright (V < 11.0) and restricted 
to the central 20" of the pickle. 

3. Proxima Cen exhibits four distinct segments with two distinct behavior modes: short 
period, low amplitude and long period, large amplitude. These variations are consistent with 
a fundamental rotation period, P = SJ^B and three darker spots. Two of the spots are either 
very small or very low-contrast. They are spaced by 180° and persist throughout our temporal 
coverage, over 24 rotations. A single more prominent spot (either large or high-contrast) formed 
in less than one rotation period, persisted through four rotations, then disappeared. A spot 
reappeared within 5° of this same longitude five rotations later. The hemisphere opposite this 
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spot brightened each time the spot formed. 

It is intruiguing that active longitudes of spot formation separated by 180° are observed 
in chromospherically active stars with close stellar companions ( Penry et al., 1995| ). If the 
photometric behavior of Proxima Cen is indicative of a synchronously rotating companion, its 
mass is less than that of Jupiter (Benedict et al., 1997a). 



4. We interpret the four distinct segments with two distinct behavior modes seen in the 
Proxima Cen photometry as an activity cycle with a period ~ 1100°'. Most of the flare activity 
occurred in the long period, high amplitude variation segment B. In the phased lightcurve three of 
the four detected flares are near the deepest minimum. 

5. The scatter in the Proxima Cen phased light curve is far larger than our photometric 
precision. This scatter could be caused by the formation and dissolution of small spots or plages 
within one rotation period. 

6. We find brightness variations five times our formal photometric precision for Barnard's 
Star. Unfortunately, the sparse coverage of the possible variation renders it a marginal detection. 
We conclude that Barnard's Star shows very weak evidence for periodicity on a timescale of 
approximately 130 days 

To confirm the spots and the inferred rotation periods will require observations of color 
changes (e.g. Vrba et al., 1988| ) and additional spectroscopic observations of lines sensitive to 



presence or absence of star spots. Extended-duration milli-magnitude V band photometry from 
the ground, while difficult ( pilliland et al., 1993 ), could probe the activity cycle periodicity of 



Proxima Cen. Future tests could include Space Interferometry Mission ( Shao, 1995D observations 
with several microarcsec astrometric precision. If spots and plages exist on these stars, we can 
expect easily detectable star position shifts as activity sites vary. Such observations will provide 
detailed maps of spot and plage location. Extended temporal monitoring will provide evolutionary 
details. 
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Appendix 1 



Appendix 1 contains the observation logs and measured average S values for Proxima Cen 
and Barnard's Star. 
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Table 1. Proxima Cen 



Parameter 


Value 


aliases 


a Cen C, G J 551, V645 Cen 


Mv 


15.45 ±0.1 


B-V 


1.94 


Sp.T. 


M5Ve 


Mprox 


O.11M0 




O.OO1L0 


-Rprox 





Reference 



Gliese &: Jahreiss, 1991 



Kirkpatrick fc McCarthy 1994 



Liebert &: Probst 1987] 



Panagi &: Mathioudakis 1993 
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Table 2. Barnard's Star 



Parameter 



Value 



Reference 



aliases GJ 699, G 140-24, LHS 57 

My 13.2 ±0.1 

B-V 1.73 

Sp.T. M4Ve 

MBarn O.17M0 

i^Barn 0.0046^0 

i?Barn O.17i?0 



Gliese & Jahreiss, 1991 



Kirkpatrick &: McCarthy 1994 



Henry fc McCarthy, 1993| 



Henry fc McCarthy, 1993 
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Table 3. Flat Field Modeling - M35 Stars 



Parameter S9_79 Sll_68 S9_47 

linear trend 

loicounts) 1.134 1.107 1.129 

±0.008 0.023 0.015 

r (counts d-i) -1.36E-05 -l.llE-05 -1.31E-05 

O.lOE-05 0.24E-05 O.lOE-05 

Variation from 180 deg flip 

A (counts) 0.0006 -0.0065 0.0039 

0.0006 0.0003 0.0004 

P (days) 365.25 365.25 365.25 



(f) (radians) 281.56 281.71 282.00 

1.03 0.17 0.14 
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Table 4. Flat Field Modeling - M Dwarf Stars 



Parameter GJ 623 and GJ 748 



Io{counts) 1.133 

±0.0132 

r (counts d-^) -1.30E-05 

0.13E-05 
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Table 5. Proxima Cen Light Curve Parameters 



Segment 


N 


lo 


/' 


A 


P 




(counts) 


(counts d~^) 


(counts) 


(days) 


(radians) 


A 


32 


3560 


0.05 


34.2 


41.8 


-1.6 






±710 


0.08 


6.0 


0.9 


2.7 


BCD 


39 


3200 


0.09 


-102 


82.5 


-0.4 






160 


0.02 


12 


0.3 


2.4 


B 


17 


3070 


0.10 


-139 


82.7 


-2.8 






580 


0.06 


9.5 


0.7 


0.7 


C 


11 


4250 


-0.02 


-34.5 


42.4 


-2.9 






200 


0.02 


8.0 


0.2 


1.3 


D 


11 


5473 


-0.14 


-114.8 


84.3 


-1.5 






645 


0.06 


16.3 


0.8 


2.2 
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Barnard's Star Reference Frame 
06 August 1994 



ID 


V 


BARN 


9.5 


30 


13.2 


31 


15.5 


32 


14.8 


33 


14.6 


34 


14.1 


35 


15.7 


36 


11.6 



'O °34 

36 



O2000 



30 
o 



.32 



35 - 1993.2 
O1990 



33 



+x 



FGS3 

\ 



+y 



Fig. 1. — The Barnard's Star field on 4 August 1994. North is at top, east to the left. The pickle 
coordinate system (x, y) is indicated. The symbol size is proportional to the relative brightness 
of each reference star (listed). The central circle (diameter ~ 3f8) is accessible by the FGS 3 
instantaneous aperture for any HST roll. Barnard's Star is labeled at three epochs. 
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Fig. 2. — S (average counts per 0.025 sec) for two faint astrometric reference stars in the Barnard's 
Star field (stars 32 and 33, Figure 0) vs. target - Sun separation in degrees. Fitting function is 
equation 0. 



- 25 - 



140 



138 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 1 1 



Lightcurve for Average of 






Barnard Reference Stars 32 & 33 






□ 

□ ° 


□ 




° J □ 

nUU^ 


□ 

■ 

• ■ 


□ 

□ 


^ □ ■ ■ *. ■ . " 


■ 


• 


■ Corrected for Zodiacal Light 






□ original 







132 - 



130 - 



9200 9400 9600 9800 
JD-2440000 (days) 



10000 



10200 



Fig. 3. — Correcting for Zodiacal Light. Average intensity (S, counts per 0.025 sec) vs. modified 
Julian Date for the average of two faint astrometric reference stars in the Barnard's Star field 
(Figure The horizontal line denotes the mean brightness of the corrected data. Scatter is far 
less in the corrected data. Uncorrected and corrected values have been flat-fielded with equation 
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Fig. 4. — Field of view of FGS 3 with position of three M35 stars. The V magnitude is given as 
part of the identification (e.g., S9_ 47 has V = 9.47). 
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Fig. 5. — Time-dependent photometric variations for three stars in M35 modeled as linear trends 
whose parameters (I and I') are given in Table ^. The residuals in the top panels are size-coded 
to show the two LTSTAB orientations (Fall = small, Spring = large) and fit with sine functions. 
Note differences in the variation amplitude (parameter A in Table |3|) . 
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Fig. 6. — Trend of normalized intensity for Gl 623 and GJ 748. The fitted line parameters are 
given m Table | 
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Fig. 7. — Amplitude of side-to-side variation vs. distance from pickle center. Targets are stars from 
M35 (location from Figure]^ and amplitude from Figure stars from the Barnard's Star reference 
frame (Figure |^), and GJ 748. Error bars along the X-axis indicate the radial range within the 
pickle for all observations of the particular target. 
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Fig. 8. — Intensities of Barnard field reference stars 34 and 36 plotted against pickle coordinates 
for all observation dates. The symbol size is proportional to the normalized intensity. The low- 
sensitivity region is clearly seen in the vicinity of (x, y) — (-25, 627). 
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Fig. 9. — Photometry of Proxima Cen. Each orbit contains 2, 3, or 4 exposures (Appendix 1.1). 
Exposures perturbed by flares are marked Fl - F4. Exposures acquired during CVZ orbits are 
labeled 'cz'. Error bars are about the size of the plotted symbols. Four segments and two distinct 
behavior modes are identified, A - D. Trend line is discussed in text. 
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Fig. 10. — Photometry of Barnard's Star. Each symbol represents the average of an exposure. 
Each orbit contains 3 exposures (Appendix 1.2). No flares were detected. Error bars are smaller 
than the plotted symbols. The time scale is identical to Figure ^. Trend line is discussed in text. 
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Fig. 11. — Periodogram from 71 normal points (average for each orbit) with flares removed, flat- 
fielded with equation ^. The periodogram has most significant peak at P ~ 83 days, with less than 
a 0.1% false-positive probability. 
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Fig. 12. — Periodogram from 35 normal points (average for each orbit), flat-fielded with equation]^. 
Periodogram has most significant peak at P = 130.4'^, with a 10% false-positive probability. 
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Fig. 13. — Direct light curve. Each symbol represents the average of from two to four exposures 
per orbit, with flares removed. Error bars are 1 — a. See Table ^, lines 1, 3, 4, and 5 for the results 
of fitting a sin wave and trend (equation ^) to each segment. 
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Fig. 14. — Phased light curves for Proxima Cen. Top: the two long-period segments (B and D, 
Figure [I^) . Sin wave fit has a period constrained to one cycle. Bottom: short-period segments (A 
and C) are phased to the longer period and show a double sin wave. The fit is constrained to have 
a period of one-half cycle. Error bars are about the size of the symbols. The observed flares are 
labeled Fl - F4. 
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Fig. 15. — Phased light curve for Proxima Cen. Bottom: the two long-period segments (B and D, 
Figure IH), phased to P = 83'?5 and fit with a sin wave having a period of one cycle. Top: the 
residuals to the sin wave fit. These show a double sin wave pattern nearly identical to the two short- 



period segments, A and C (Figures and |14|), suggesting that the low-amplitude, short-period 
signature persists during segments B and D. 
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Fig. 16. — Per-orbit average direct light curve. Each symbol represents the average of three 
exposures. Error bars are 1 — a. Bottom: fit with a trend line. Top: residuals to trend line, fit 
with a constant-amplitude sin wave with constrained P = 130'?4. 
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Fig. 17. — Trend-corrected light curve for Barnard's Star, phased to P = 130^4. The full-amphtude 
of the variation is ~ 0.01 magnitude. Error bars are 1 — a. 



